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SECTION  I 
INTRODUCTION 


Certain  experiments  and  tests  require  the  measurement  of 
out-of-plane  motions  of  target  surfaces.  One  way  these  measure¬ 
ments  can  be  made  is  with  use  of  special  devices  which  produce 
moire  fringe  patterns  on  the  surface  of  the  targets.  A  descrip¬ 
tion  of  one  of  these  devices  has  been  provided  in  a  previous 
document.  This  document  provides  detailed  descriptions  of 
procedures  related  to  the  use  of  that  device  and  the  analysis 
of  the  resulting  fringe  patterns. 

Although  the  device  has  been  used  primarily  to  determine 

the  response  of  impulsively  loaded  targets  of  fairly  large 
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area  (1500  to  2000  cm  ) ,  information  provided  in  this  document 
also  applies  to  analysis  of  data  from  surfaces  of  smaller  area 
and  surfaces  subjected  to  slowly  applied  or  static  loads.  In 
all  cases,  deflection  data  are  obtained  after  analysis  of 
photographs  of  individual  fringe  patterns. 

A  computer  program  named  PANDIS  is  used  to  assist  in  the 
reduction  of  fringe  pattern  data  to  displacement  data.  This 
program,  written  in  the  BASIC  language,  is  designed  for  use  from 
a  remote  time-shared  computer  terminal.  Section  II  of  this 
document  provides  a  brief  description  and  instructions  for  use 
of  the  device.  Sections  III,  IV,  and  V  provide  detailed  data 
reduction  instructions,  a  description  of  program  PANDIS  (in¬ 
cluding  a  listing  of  program  statements) ,  and  a  sample  problem 
and  the  resulting  printout,  respectively. 


(1)  Piekutowski,  A.  J.,  A  Device  to  Determine  the  Out-of-Plane 
Displacement  of  a  Surface  Using  a  Moire  Fringe  Technique, 
AFWAL-TR-81-3005 ,  Air  Force  Wright  Aeronautical  Laboratories, 
Wright-Patterson  Air  Force  Base,  OH,  May  1981. 


SECTION  II 
USE  OF  MOIRE  DEVICE 

As  constructed,  the  moire  device,  shown  schematically  in 
Figure  1,  simply  generates  a  system  of  light  and  dark  bands 
which  are  fixed  in  space  relative  to  the  moire  device.  The 
light  and  dark  bands  are  apparent  to  the  viewer  only  when  the 
system  of  bands  is  interrupted  by  a  target  surface.  A  critical 
orientation  of  the  target  surface  with  respect  to  the  moire 
device  is  not  required  since  the  initial  shape  of  the  target 
surface  can  always  be  reconstructed  after  analysis  of  the  fringe 
patterns . 

An  illustration  of  the  technique  used  to  produce  the  moire 
patterns  is  presented  in  Figure  2.  As  shown  in  this  figure, 
the  optical  system  on  the  right  (Projection  System)  projects 
an  image  of  the  Ronchi  ruling  in  the  system  onto  the  target 
surface.  The  lens  in  the  left  optical  system  (Viewing  System) 
focuses  an  image  of  the  target  and  projected  ruling  onto  the 
Ronchi  ruling  in  the  viewing  system.  A  moire  pattern  is  pro¬ 
duced  where  the  image  of  the  projected  ruling  is  superimposed 
on  the  viewing  system  Ronchi  ruling.  A  camera,  with  its  lens 
focused  on  the  composite  image,  is  used  to  photographically 
record  the  moire  pattern. 

Analysis  of  raw  data  is  facilitated  by  exercising 
reasonable  care  during  setup  of  the  moire  device.  Attention 
is  drawn  to  the  relationship  of  the  shaded  regions  adjacent 
to  the  rays  (boundaries  between  clear  and  opaque  areas  of 
Ronchi  ruling)  shown  in  Figure  2.  The  regions  shown  are 
adjacent  to  the  optical  axes  of  the  projection  and  viewing 
systems.  When  the  rays  associated  with  these  regions  and  the 
target  surface  intersect  simultaneously ,  a  dark  band  will 
appear  at  the  target  surface.  This  particular  dark  band  is 
identified  as  fringe  zero.  Identification  of  the  rays  which 
produce  fringe  zero  is  facilitated  by  heavy  locator  lines 


BOTH  SYSTEMS  NOT 


NEED  TO  BE  AT  SAME  ANGLE 


TO  THE  TARGET  SURFACE 


placed  on  the  Ronchi  rulings.  Relative  positions  of  the  locator 
lines  in  the  various  dark  bands  of  the  moire  pattern  are  used  to 
establish  the  direction  which  the  moire  device  must  be  moved  to 
produce  the  desired  target-device  relationship.  The  positions 
of  the  locator  lines  are  shown  for  a  number  of  fringe  surfaces 
in  Figure  3.  The  coordinate  system  and  fringe  numbering  con¬ 
vention  which  was  adopted  for  use  during  subsequent  data  reduc¬ 
tion  procedures  are  also  shown  in  this  figure. 


PROJECTION  OF  OPTICAL  AXIS 


Figure  3.  Arbitrary  System  of  Identifying  Fringes  and  Coordi¬ 
nate  System  Used  in  Analysis  of  Fringe  Patterns.  In 
this  figure,  the  reference  plane  is  parallel  to  the 
X-D  plane.  Note  the  position  of  locator  lines  in 
each  fringe. 


Information  presented  in  Figure  4  is  used  to  illustrate 
relationships  between  the  various  device  parameters  which  must 
be  determined  for  correct  interpretation  of  the  resultant 
fringe  patterns.  However,  one  important  difference  between 
the  device  and  the  illustration  in  Figure  4  exists.  For  pur¬ 
poses  of  simplifying  presentation,  the  Ronchi  ruling  has  been 
shown  between  the  appropriate  lens  and  the  target  rather  than 
behind  the  lens,  e.g.,  Figures  1  and  2.  Positioning  the  Ronchi 
ruling  as  shown  in  Figure  4  merely  facilitates  a  graphical 
presentation  of  the  rays  used  to  determine  the  location  of  the 
dark  bands.  For  consistency  of  illustration,  the  lens  on  the 
right  side  of  Figure  4  corresponds  to  the  viewing  lens  of 
Figure  2  (on  the  left) . 

Selection  of  nominal  characteristics  of  various  optical 
components  used  in  the  device  is  made  during  design  of  the 
device.  Determination  of  actual  characteristics  of  these 
components  is  made  during  alignment  and  assembly  of  the  device. 
The  six  parametric  constants  determined  during  assembly  of  the 
device--aT ,  a_,  RT ,  R0,  PT ,  and  P  --are  shown  in  Figure  4. 

li  K  Li  K  Li  K 

These  constants  are  used  to  uniquely  specify  the  locations  of 
the  various  dark  bands  produced  by  the  moire  device.  Following 
is  a  brief  description  of  each  of  these  parameters.  Subscripts 
R  and  L  refer  to  the  right  and  left  side  of  Figure  4. 


Parameter 


Description 

Angle  optical  axis  makes  with  X-axis 

Normal  distance  from  center  of  lens  to 
X-axis  or  target  surface  (or  tangent  to 
target  surface  in  case  of  curved  targets) 
In  actual  practice  this  value  is  calcu¬ 
lated  using  the  relationship  shown  in 
Figure  4  and  the  appropriate  values  of  a 
and  the  distance,  d,  which  was  measured 
during  alignment  of  the  individual  heads. 


Parameter 


Description 

P  A  pseudopitch  used  when  determining,  in 

slope-intercept  form,  the  equations  of 
lines  representing  the  various  rays.  This 
dimensionless  value  is  determined,  as 
shown  in  Figure  4,  using  the  Ronchi  ruling 
pitch  (furnished  by  Ronchi  ruling  manufac¬ 
turer)  ,  the  system  magnification  (measured 
during  the  alignment  of  the  individual 
heads) ,  and  the  previously  mentioned 
distance,  d. 

Target  surfaces  shown  in  the  four  figures  just  presented 
are  represented  as  planar  surfaces  which  are  tangent  to  fringe 
zero.  This  orientation  is  not  required  for  satisfactory  opera¬ 
tion  of  the  device,  and  practically  speaking,  is  rarely  achieved. 
When  the  device  is  used  with  convex  or  concave  surfaces,  the 
preferred  orientation  of  the  device  has  the  displacement  axis 
of  Figure  3  coincide  with  a  line  normal  to  a  plane  tangent  to 
the  target  surface  at  the  point  where  intersection  of  the  X-axis 
interruptions  in  the  locator  lines  occurs.  The  relationship 
of  the  moire  device  displacement  axis  (fixed  by  optical  compo¬ 
nents  in  the  moire  device)  and  a  line  normal  to  the  target  at 
the  point  of  tangency  just  described  is  the  preferred  orienta¬ 
tion  of  the  device  and  target  for  all  applications  since  it 
facilitates  identification  of  fringe  zero.  When  the  device  is 
mounted  independent  of  the  target,  the  user  must  carefully 
position  the  device  for  each  test  to  insure  this  orientation 
is  achieved  as  closely  as  possible.  In  other  applications, 
e.g.,  measurement  of  aircraft  transparency  deflections,  the 
moire  device  may  be  installed  in  mounts  which  may  or  may  not 
provide  for  adjustment  of  the  relationship  of  the  device  to 
the  target  surface.  In  these  cases,  the  general  relationship 
of  the  target  and  the  device  is  fixed  during  design  of  the 
device  mounting  hardware . 

Most  uses  of  the  moire  device  require  some  preparation 
of  the  target  surface.  Typical  preparations  consist  of  the 
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application  of  a  reflective  coating,  usually  white  paint,  and 
the  installation  of  fiducial  marks  on  the  reflective  screen. 

The  fiducial  marks  are  used  to  provide  data  for  film  magnifica¬ 
tion  determinations.  A  suggested  fiducial  line  layout  is 
presented  in  Figure  5.  Use  of  X*s  rather  than  short  line  seg¬ 
ments  as  reference  points  on  the  horizontal  axis  eliminates 
the  likelihood  of  one  or  more  of  the  line  segments  being 
obscured  by  the  lines  in  the  projected  or  viewed  Ronchi  ruling. 


DIRECTION  OF  LINES 
IN  RONCHI  RULING 


2  INCH  SPACING  TYPICAL 


IDEAL  CENTER  OF  MOIRE  DEVICE 

V 

,  V.  V,  V 
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2 INCH  SPACING  TYPICAL 


COINCIDENT  WITH  X-AXIS^. 
OF  MOIRE  DEVICE 


NUMBER  ANO  LOCATION  OF  FlOUClAL 
MARKS  CAN  BE  VARIED  TO  SUIT 
TEST  CONDITIONS 


Figure  5.  Suggested  Fiducial  Pattern  for  Use  with  Moire  Device 


SECTION  III 
DATA  REDUCTION 

Data  used  to  generate  deflection  information  are  obtained 
from  various  photographs  taken  during  the  experiments  or  tests. 

A  single  photograph  may  be  analyzed  in  the  case  of  a  statically 
loaded  target.  When  the  target  is  impulsively  loaded,  data  are 
taken  from  photographs  of  selected  frames  of  a  high-speed  film 
of  the  transient  fringe  patterns.  Measurements  of  the  locations 
of  timing  marks  placed  on  the  high-speed  films  are  used  to 
determine  the  framing  rate  of  the  camera  and  the  time  interval 
between  successive  fringe  patterns  for  these  dynamic  events. 

In  this  section,  detailed  instructions  for  making  a  variety  of 
measurements  are  provided.  The  measurements  are  described  and 
presented  in  the  order  in  which  they  are  entered  into  the  com¬ 
puter  program.  A  listing  and  description  of  sample  data  items 
is  presented  in  Section  V.  The  first  data  entries  consist  of 
specific  moire  device  parameters  and  information  unique  to  the 
test.  An  explanation  of  each  of  these  data  items  is  given 
in  Section  V. 

Although  the  original  negative  film  may  be  used  for  data 
reduction  purposes,  analysis  of  the  moire  patterns  is  greatly 
facilitated  by  using  enlarged  prints  of  the  individual  frames. 

In  the  case  of  multiple  prints,  it  is  extremely  important  that 
all  frames  taken  from  a  data  film  be  printed  at  the  identical 
magnification.  Data  used  to  obtain  film  magnification  infor¬ 
mation  are  taken  from  pre-event  photographs.  Post  event 
photographs  are  used  to  obtain  measurements  of  fringe  location 
and,  ultimately,  target  displacement.  Since  the  fiducial  lines 
are  applied  to  a  surface  which  is  in  motion  during  the  event, 
film  magnification  information  taken  from  post  event  frames 
should  be  regarded  as  suspect. 

Because  the  axis  of  the  viewing  system  is  at  an  angle  to 
the  target  surface,  distortion  of  the  image  produced  at  the 
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viewing  system  Ronchi  ruling  occurs.  Consequently,  magnifica¬ 
tion  of  the  target  image  varies  continuously  in  the  horizontal 
(X)  and  vertical  (Y)  directions.  Since  out-of-plane  displace¬ 
ment  data  are  derived  from  measurements  taken  from  the  film, 
conversion  of  measured  dimensions  (film)  to  true  dimensions 
(target)  requires  that  the  variation  in  film  magnification  be 
accounted  for  properly. 

Film  magnification  along  the  X-axis  is  determined  using 
measurements  taken  from  one  or  more  of  the  frames  immediately 
preceding  the  event.  (Only  one  set  of  averaged  values  is 
entered  into  the  program,  however.)  Profiles  other  than  those 
along  the  X-axis  are  usually  taken  from  sections  parallel  to 
the  X-axis  and  through  the  short  horizontal  lines  on  the  two 
long  vertical  lines  of  the  fiducial  system  (Figure  5) .  Mag¬ 
nification  in  the  Y  direction  is  thus  treated  graphically  and 
no  further  corrections  are  required.  An  illustration  of  a 
modified  target-viewing  head  relationship  is  presented  in 
Figure  6.  In  this  figure,  film  normally  in  the  camera  is  shown 
in  the  plane  of  the  viewing  head  Ronchi  ruling.  A  view  of  the 
undistorted  target  surface  and  fiducial  system  is  shown  in 
the  upper  right  of  this  figure.  A  view  of  the  distorted  image 
produced  on  the  film  is  shown  at  the  lower  right  of  Figure  6. 
Also  shown  in  this  figure  is  the  position  of  the  interruption 
in  the  locator  line  with  respect  to  the  X-axis. 

As  shown  in  the  lower  right  portion  of  Figure  6,  the 
locator  line  (ray  zero)  does  not  coincide  with  the  center  of 
the  fiducial  pattern  and  an  offset  of  unknown  magnitude  results. 
The  magnitude  of  this  offset  and  actual  film  magnification 
data  are  determined  graphically  using  measured  values  taken  from 
the  pre-event  photographs.  All  measurements  are  made  with 
respect  to  the  locator  line  of  the  viewing  head.  This  line 
will  always  appear  as  a  straight  line  in  all  photographs.  Raw 
film  magnification  data  consist  of  a  table  of  measured  distances 
to  known  points  (X's  on  the  horizontal  fiducial  line)  and  the 
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Figure  6.  iiius cracion  ot  Technique  wnich  Produces  a  Distorted 
Image  of  Target  Surface  and  Fiducial  Pattern.  Refer¬ 
ence  to  Figures  3  and  4  will  clarify  the  apparent 
reversal  of  right  and  left  in  the  distorted  image  of 
the  target  surface  and  fiducial  pattern. 

corresponding  true  distance  to  these  points.  Measured  distances 
to  the  right  of  the  locator  line  have  positive  values  and  mea¬ 
sured  values  to  the  left  of  the  locator  line  have  negative 
values.  When  these  data  are  plotted  as  shown  in  Figure  7,  the 
magnitude  of  the  offset  between  the  center  of  the  moire  device 
and  the  center  of  the  fiducial  pattern  corresponds  to  the  true 
distance  between  the  parallel  lines  representing  the  center  of 
the  moire  device  and  the  center  of  the  fiducial  system. 

The  value  of  offset  thus  obtained  is  used  to  "correct" 
the  true  values  and  generate  a  table  of  "corrected"  film  mag¬ 
nification  data  which  are  entered  into  program  PANDIS.  For 
example,  the  true  position  of  the  second  X  to  the  right  of 
center  of  the  fiducial  pattern  is  +4  inches.  This  true  value 
is  corrected  by  subtracting  the  algebraic  value  of  the  offset 
from  the  true  position  of  the  fiducial  mark,  i.e.,  +4  -  (+1.2)  = 
+2.8.  All  other  true  values  are  corrected  in  a  like  manner. 
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LOCATION  OF  X's  ON  HORIZONTAL  FIOUCIAL  LINE  (INCHES) 

Figure  7.  Illustration  of  Graphical  Method  Used  to  Determine 
Offset  and  "Corrected"  Film  Magnification  Data. 

The  measured  distance  from  the  locator  line  to  this  X  (and  all 
other  X’s  in  the  photograph)  remains  unchanged. 

Processing  of  raw  film  magnification  data  in  PANDIS  yields 
the  coefficients  of  a  second  order  polynomial.  Subsequent  con¬ 
version  of  measured  distances  taken  from  data  frames  to  true 
distances  is  accomplished  by  evaluating  the  polynomial  using 
the  measured  values  as  the  independent  variable. 

The  framing  rate  of  the  high-speed  camera,  in  pictures  per 
second  (pps) ,  is  determined  for  that  portion  of  the  high-speed 
film  in  which  response  of  the  target  is  recorded.  Framing  rate 
is  determined  from  measurements  illustrated  in  Figure  8.  In 
this  figure,  timing  mark  zero  was  arbitrarily  selected  as  the 
sixth  mark  before  the  frame  in  which  response  of  the  target  was 
first  observed.  In  practice  any  of  the  several  marks  appearing 
before  this  frame  may  be  selected  as  timing  mark  zero.  All 
subsequent  measurements  are  made  with  respect  to  timing  mark 
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FIRST  FRAME  SHOWING  EVENT 


TIMING  ^ 

MRK  0  .TIMING  MARKS  AT  1  msec. 

|  INTERVALS  \ 

LAST  FRAME  READ\ 

rl 

\7T\ 

i  •  « 

1 

)  _  •  •  \ 

j  »  »  #  \ 

um*o 

*U(2  )*l 

- - U(3) - - 

- - U(4 

— U(0) 

J 

[ 

1 - * 

1  ENLARGED  VIEW  OF 

J  SPROCKET  HOLE 

J  (NOTE  SIDE  OF  SPROCKET 
]  HOLE  WHICH  IS  USED) 

Figure  8.  Illustration  of  a  Portion  of  High-Speed  Film  Showing 
Measurement  Points  for  Film  Speed  Determination. 

zero.  In  program  PANDIS,  the  time  to  successive  marks,  R(  )  , 
and  the  corresponding  measured  distances,  U(  ),  are  entered 
as  data  pairs  and  used  in  conjunction  with  certain  physical 
characteristics  of  16  mm  film  to  determine  the  coefficients 
of  a  polynomial  which  yields  the  instantaneous  framing  rate 
of  the  camera  as  a  function  of  frame  number. 

The  frame  of  film  immediately  before  the  frame  in  which 
response  of  the  target  is  observed  is  identified  as  frame  zero. 
Frames  exposed  after  frame  zero  are  numbered  sequentially  in 
ascending  order.  The  first  measurement  to  be  taken,  U(0),  is 
used  to  determine  the  location  oi;  frame  zero  with  respect  to 
timing  mark  zero.  The  other  measured  distances,  U(l),  U( 2), 
U(3),  etc.,  are  used  to  determine  the  actual  framing  rate  of 
the  camera.  Note  that  the  distance  to  timing  mark  zero,  U(l), 
is  zero  and  entered  accordingly  in  statement  100  of  the  Sample 
Data  (Section  V) . 


Since  the  system  of  dark  bands  generated  by  the  moire 
device  is  unique  for  the  set  of  system  parameters,  displacement 
of  a  target  at  any  point  can  be  determined  if  the  X-coordinate 
and  fringe  number  of  the  point  are  specified.  The  X-coordinate 
and  fringe  number  are  obtained  from  analysis  of  a  photograph 
of  the  fringe  pattern.  A  simulated  photograph  of  an  ideal 
target  and  fringe  pattern  is  presented  in  Figure  9.  The  intended 
use  of  the  out-of-plane  displacement  data  determines  the  number 
of  measurements  taken  from  the  photograph.  If  the  deflected 
shape  along  one  section  is  required  (usually  along  the  horizon¬ 
tal  fiducial  line) ,  measurements  are  taken  where  the  dark  bands 
intersect  the  fiducial  line.  If  a  map  of  the  deflected  surface 
is  desired,  additional  profiles  are  obtained  along  sections 
parallel  to  the  horizontal  fiducial  line.  Examples  of  both 
types  of  measurements  are  illustrated  in  Figure  9. 

Analysis  of  photographs  begins  with  identification  of 
the  various  dark  bands  in  the  moire  pattern.  Use  of  the  locator 


Figure  9.  Simulated  Photograph  of  Ideal  Moire  Pattern  Showing 
Various  Measurement  Locations. 


lines  aids  in  the  identification  of  fringe  zero.  Study  of 
Figures  3  and  4  will  assist  in  the  identification  of  all  other 
dark  bands  in  the  pattern.  If  deflected  surface  data  are  re¬ 
quired,  lines  are  drawn  on  the  photograph  through  the  short 
horizontal  lines  which  are  part  of  the  vertical  fiducial  lines. 
Finally,  measurement  of  the  various  distances  is  made.  As  shown, 
these  distances  are  measured  from  the  center  of  the  dark  band 
to  the  locator  line.  Note  that  all  measurements  are  made  per¬ 
pendicular  to  the  locator  line  and  that  the  correct  algebraic 
sign  of  the  measured  value  must  be  observed. 

Data  for  the  measured  points  along  each  section  line  are 
entered  as  fringe  number-measured  value  pairs  beginning  with 
the  pair  at  the  left  of  the  photograph  and  ending  with  the  pair 
at  the  right  of  the  photograph.  When  data  for  a  number  of 
sections  from  a  photograph  are  to  be  entered,  data  from  the 
uppermost  section  are  entered  first;  order  of  entry  proceeds 
successively  to  the  lower  sections. 

Procedures  used  in  the  collection  of  fringe  number- 
measured  value  data  are  identical  for  the  pre-event  and  post 
event  photographs  with  one  exception.  An  additional  measurement 
is  taken  from  the  first  post  event  frame  or  the  pre-event  frame 
if  the  target  is  suitably  identified.  This  measurement,  XO  in 
Figure  9,  is  the  distance  from  the  locator  line  to  the  point 
of  the  event  or  the  center  of  the  area  of  interest.  The  value 
of  XO  is  entered  into  PANDIS  prior  to  entry  of  the  pre-event 
and  post  event  data. 

Examination  and  study  of  the  sample  data  in  Section  V 
should  clarify  and  further  the  program  user's  understanding  of 
data  entry  procedures  used  with  PANDIS. 


SECTION  IV 
PROGRAM  DESCRIPTION 

1.  PROGRAM  USAGE 

Program  PANDIS  consists  of  a  main  program  and  two  sub¬ 
routines,  LSQ  and  COR.  The  main  program  directs  reading  of 
input  variables,  raw  data,  order  of  processing  raw  data,  and 
printing  of  results.  Subroutine  LSQ  is  called  upon  several 
times  during  execution  of  the  main  program  to  fit  data  sets  to 
least-squares  polynomials.  Subroutine  COR  is  called  after  each 
raw  data  pair  has  been  read  and  processed  to  yield  a  displace¬ 
ment  value  unique  to  the  fringe  number  and  X-coordinate  of  the 
data  pair.  In  subroutine  COR,  the  initial  value  of  displacement 
is  examined  and  corrected,  if  necessary,  for  magnification 
errors  which  are  introduced  when  the  target  surface  (hence, 
measurement  plane)  has  been  displaced  a  sufficient  distance  from 
the  initial  reference  (pre-event)  plane.  A  simplified  flow 
chart  for  the  main  program  is  presented  in  Figure  10. 

Input  variables  and  test  peculiar  data  are  read  at  the 
beginning  of  the  main  program.  Certain  of  these  data  are  used 
to  characterize,  in  tabular  form,  the  location  (in  displacement, 
D,  and  X-coordinates)  of  the  fringes  which  were  observed  in  the 
test  film.  In  the  next  portion  of  the  program,  data  used  to 
determine  film  magnification  and  speed  are  processed  to  yield 
least-squares  polynomials  which  describe  the  variations,  with 
distance  and  time,  respectively,  of  these  film  parameters.  Work 
on  actual  fringe  pattern  data  begins  with  data  taken  from  the 
pre-event  frame.  Individual  sets  of  X-D  data  taken  from  this 
frame  are  fit  to  least-squares  polynomials  which  attempt  to  des¬ 
cribe  the  shape  of  the  undeflected  target  surface  along  a  partic* 
ular  measurement  section  line.  At  the  completion  of  printout  of 
the  results  of  the  first  fit,  the  user  is  asked  if  the  next 
higher  order  of  fit  should  be  made.  A  positive  response 
results  in  another  pass  through  subroutine  LSQ.  Printout  of 
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Figure  10.  Flowchart  for  Program  PANDIS.  Initial  values  for 

range  of  program  statement  numbers  in  which  specific 
operations  are  performed  are  shown  at  left  of  flow¬ 
chart  symbols. 
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the  results  of  this  fit  is  again  followed  by  the  opportunity  to 
request  the  next  higher  order  of  fit.  After  various  orders  of 
fit  have  been  examined,  the  user  is  asked  to  select  that  order 
of  polynomial  which  best  describes  the  pre-event  panel  shape. 
Upon  receipt  of  a  reply,  control  of  processing  of  data  from 
post  event  frames  and  the  printing  of  results  is  returned  to 
the  main  program. 


2.  DESCRIPTION  OF  VARIABLES 

Following  is  a  listing  (and  brief  description)  of  the 
simple  variables,  arrays,  and  string  variables  used  in 
program  PANDIS.  Where  appropriate,  units  and  applicable 
reference  figures  are  included  in  this  listing. 

Simple  Variables 

A1  -  Angle,  a^,  which  optical  axis  of  R.H.  optical  system 
makes  with  X-axis,  radians  (Figure  4) . 

A2  -  Angle,  otL/  which  optical  axis  of  L.H.  optical  system 
makes  with  X-axis,  radians  (Figure  4) . 

A3  -  Angles  which  rays  in  R.H.  optical  system  make  with 
X-axis,  radians. 

A4  -  Angles  which  rays  in  L.H.  optical  system  make  with 
X-axis,  radians. 

A5  -  Angle,  a,  used  when  "correcting"  displacements  for 
magnification  errors,  radians  (Figure  11)  . 

A6  -  Angle,  3,  used  when  "correcting"  displacements  for 
magnification  errors,  radians  (Figure  11) . 

E  -  Standard  error  of  estimate  of  least-squares  fit. 

FO  -  Time  after  event,  evaluated  for  each  frame  number, 
msec . 

FI  -  Largest  positive  fringe  number  used  in  evaluation. 

F2  -  Largest  negative  fringe  number  used  in  evaluation. 

F3  -  Average  framing  rate  of  camera  over  an  arbitrary 

1  msec  interval,  IQ,  pps . 
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Average  framing  rate  of  camera  over  the  arbitrary 
1  msec  interval ,  I  +  1 ,  pps . 

Aq  term  of  polynomial  describing  film  speed. 

Ai  term  of  polynomial  describing  film  speed. 

A2  term  of  polynomial  describing  film  speed. 

General  operator. 

Operator  to  sequence  through  all  frames  of  film  after 
pre-event  frame. 

Operator  to  sequence  through  all  sections  of  a  frame. 
General  operator. 

General  operator. 

General  operator. 

General  operator. 

Order  of  fit  desired  in  least-squares  subroutine. 

Ml  +  1. 

Aq  term  of  polynomial  describing  film  magnification, 
inches. 

Ai  term  of  polynomial  describing  film  magnification. 

A2  term  of  polynomial  describing  film  magnification, 
inches"! . 

General  operator. 

Number  of  frames  of  film  to  be  evaluated  (including 
pre-event  frame. 

N9  +  1. 

Frame  number. 

Number  of  sections  in  a  frame. 

Value  used  to  signal  that  insufficient  data  are 
available  for  a  least-squares  fit. 

Number  of  rays  on  either  side  of  optical  axis 
(Figure  4) . 


PO  -  Proportionality  factor  used  during  interpolation  in 
X-D  Tables. 

PI  -  Number  of  data  pairs  for  a  particular  use  of  least- 
squares  subroutine  or  other  functions  requiring 
knowledge  of  number  of  data  pairs  for  a  particular 
operation . 

P7  -  Page  counter. 

P9  -  Pi  (3.14159) . 

R1  -  Normal  distance  from  R.H.  optical  system  lens  to 
X-axis,  mm  (Figure  4). 

R2  -  Normal  distance  from  L.H.  optical  system  lens  to 
X-axis,  mm  (Figure  4). 

R3  -  Constant  used  when  "correcting”  displacements  for 
magnification  errors,  mm. 

R8  -  Pseudopitch  of  R.H.  optical  system  (Figure  4). 

R9  -  Pseudopitch  of  L.H.  optical  system  (Figure  4). 

S  -  Summation  in  least-squares  subroutine. 

SO  -  Calculated  displacement  (w/r  to  Moire  Device)  using 
pre-impact  data  fit,  mm. 

S2  -  Ax  in  "correction"  of  displacements  for  magnification 
errors,  mm. 

56  -  Temporary  storage  of  X  values,  mm. 

57  -  Temporary  storage  of  D  values,  mm. 

XO  -  Point  of  impact  w/r  to  centerline  of  Moire  Device, 
inches  and  mm  (Figure  9) . 

YO  -  Counter  in  correction  subroutine. 

Y2  -  Summation  in  least-squares  subroutine. 

Y8  -  Projectile  weight,  variable. 

Y9  -  Impact  velocity,  variable. 

Z  -  Absolute  value  of  (Y2-S) . 
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)  -  Coefficients  obtained  in  least-squares  subroutine. 

-  Intercepts  of  rays  in  right  optical  system,  mm. 

-  Intercepts  of  rays  in  left  optical  system,  mm. 


Arrays 
A  (  , 

B  (  ) 

C  (  ) 

D  (  ,  ) 

E  (  ,  ) 

F  (  ) 

G  (  ) 

H  (  ,  ) 

M  (  ) 

N  (  ) 

0  (  ,  ) 

P  (  ) 

Q  (  ) 

R  (  ) 

T  (  ,  ) 

U  (  ) 

V  (  ,  ) 


Displacement  values  in  table  of  fringe 
coordinates,  mm. 

Specific  displacement  values  from  pre-event 
frames,  mm. 

Counter  number  for  X  and  D  values  tabulated 
for  each  fringe. 

Storage  for  various  Xi  summations  in  least- 
squares  subroutine. 

Specific  fringe  numbers  from  pre-event  frames. 

Slopes  of  rays  in  right  optical  system. 

Slopes  of  rays  in  left  optical  system. 

Storage  of  coefficients  of  polynomials  des¬ 
cribing  each  section  in  pre-event  frame. 

0(9,  )  =  PI  and  0(10,  )  =  Ml. 

X-values  in  least-squares  subroutine,  inches 
or  mm;  measured  values  for  film  magnification 
data,  inches;  and  measured  values  for  post 
event  frames,  inches. 

Y-values  in  least-squares  subroutine,  inches 
or  mm;  true  values  for  film  magnification 
data,  inches;  and  fringe  numbers  for  post 
event  frames. 

Time  values  for  film  speed  determination, 
msec  and  measured  value  for  pre-event  data, 
inches . 

Various  Exi^i  in  least-squares  subroutine. 

Measured  distances  for  film  speed  determination 
inches  (Figure  8) ,  and  fringe  number  for  pre¬ 
event  data. 

Inverse  of  Y (  ,  )  in  least-squares  subroutine. 
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W  (  ,  )  -  Specific  X-values  from  pre-event  frame,  mm. 

X  (  ,  )  -  X-values  in  table  of  fringe  coordinates,  mm. 

Y  (  ,  )  -  Various  Ex^m  in  least-squares  subroutine. 

String  Variables 

D$  -  Date. 

S$  -  Shot  or  test  identification  code. 

X$  -  Query  for  next  order  of  fit  in  least-squares 

fit  for  pre-event  data. 

P$  (  )  -  P$(0)  -  Target  description. 

P$(l)  to  P$ ( 9)  -  ST,  ND,  RD,  TH,  TH,  TH,  TH, 
TH,  TH  suffixes  added  to  order  of  fit  to  yield 
1st,  2nd,  etc. 

P$(10)  -  Projectile  description. 

P$(ll)  to  P$(20)  -  Section  location 
designators . 

3.  PROGRAM  LISTING 

Following  is  a  listing  of  the  BASIC  statements  which 
comprise  program  PANDIS .  The  program,  as  listed,  was  run  on 
a  UNIVAC  90/80-3  computer  using  the  sample  data  included  in 
the  next  section  to  produce  the  printout  presented  at  the  end 
of  the  next  section.  The  listing  has  been  separated  into  a 
number  of  statement  groups.  In  general,  each  group  corresponds 
to  one  of  the  operations  shown  on  the  flowchart  in  Figure  10. 
Where  appropriate,  a  short  explanation  of  the  function  of  each 
group  of  statements  is  presented  prior  to  the  listing  of 
statements . 

0  REM  STATEMENT  NUMBERS  1  THROUGH  999  ARE  AVAILABLE  FOR  DATA 

See  Sample  Data  (Section  V)  for  format  requirements. 


The  following  statements  allocate  storage  for  the  variously 
named  arrays. 


1  COO  II I  h  A  <  1 0  *  1 )  *  R  <  200  >  *  C  (  200  )  *  D  <  1 00 * 30  >  *  F  <  30  *  1  0  > 

1010  M  TM  f  <?0(>)  *  G  < 1 0 ) *  H  <  30  *  1 0 ) *  M  <  200 )  *N(  200)  fOdO^lO) 

i  030  n  i  m  p  *  so  >  *  r *  <  20  >  * 0  <  so ) *  r  < so  >  *  r  <  1 0  *  1  > 

1030  DIri  n  <  SO  >  10*  1  o  >  *  U  i  30*10)  *X  ( 100*30)  *Y< 10* 10) 

7000  RF  M  PARAMETRIC  COWS  FANT:;  FOR  rIOIRE  AFT1  AC  AH  IS 

20 1 0  READ  R 1  *  R2 *  N9 *  R8 *  R9 

2020  R LA U  A 1  *  A2  *  F 1  * F2 *  N3 

2030  P9-3, 14159 

2040  N1  ~  N<:M  L 

2050  FOR  1-1  TO  9 

2060  Rf  All  ft  <t) 

2070  NEXT  1 

2000  R3:=R1  /  TAN  <  A  I  > 

2090  FOR  l  ••■■■■  1  1  10  N31  10 

2100  READ  P*<I) 

3110  NEXT  1 

2200  RF  M  -  READ  GCNFRAI  INFORiiAT  1 ON  REGARDING  SHOT 
2210  READ  S*»i:i*»P*<0> 

7270  READ  P$ ( 10 ) * YO * Y9*N0 

In  the  next  operation#  the  general  form  of  equation  used 
to  determine  the  slopes  of  the  various  rays  is: 


m  =  tan  0 

where  m  is  the  slope  of  the  ray  and  0  is  the  angle  which  the  ray 
makes  with  the  X  axis.  The  D-intercept#  b#  of  each  ray  is 
computed  using  the  following  relationship  and  values  illustrated 
in  Figure  4. 


b  =  R  -  (tshO  m 


In  the  piogram,  care  is  taken  to  insure  that  the  proper  algebraic 
signs  are  observed  during  computation  of  the  various  values. 


2400  RFM  ■  DETERMINE  SLOPES  AND  INTERCEPTS  OF  ALL  RAYS 

2410  FOR  I*-N9  TO  N9 

3420  J»I+N1 

7430  A3»ATN<I*R8> 

2440  A4-ATN ( I *R9 ) 

2450  M <  J )  “ T AN <  A 1  •  A3  > 

2460  B<J)~R1  <  ( R .1  /  F AN f  A 1  )  ) #M  v  J )  ) 

2470  N<  J)a-TAN(A?-FA4) 

7480  C ( J ) =R2+ < < R2/TAN < A2 ) ) <  J ) > 

2490  NF:!XT  I 
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Before  the  coordinates  of  the  intersections  of  rays  are 
tabulated,  the  fringe  number  of  the  intersection  is  checked  to 
determine  whether  or  not  it  falls  within  the  range  of  fringe 
numbers  that  appeared  on  the  data  photographs.  Intersections  for 
fringe  numbers  which  lie  outside  the  range  are  not  computed. 
Computation  of  the  fringe  numbers  and  X  and  D  coordinates  of 
the  various  intersections  is  made  using  relationships  given  in 
Figure  4. 
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1  <  K  .»  F  <  f\  ) 

1  1 

27  A  0 
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Rl  M 
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U'-IO 

1  isr 

US  J  NO 

VO  1  0 

$07o 

1  1  S  1 

OS  1  NO 

9  c*  70 

30  30 

1.  J M 

1  IS  1  NO 

9020  ? 

S‘J> 

3040 

i  i  r.  r 

USING 

1  >030  * 

It* 

3 0  3  0 

1  JM 

UMNO 

9040  * 

J  •  F  ».  • 

30.  >0 

F.  1  S  ) 

UMNO 

’  *000  * 

Ri(  1 

30  7() 

1  IS  J 

OM  NO 

V  0  7,0  * 

ii\ 

3000 

i.  r  i 

UMNO 

Vi.i  70  «. 

i 9 

3090 

LIST 

Ob  I  NO 

9000, 

NO  | 

31  00 

II  ST 

OS  IMG 

9090 

31 10 

LIST 

OS  TNG 

°  1 00 

3:1.30 

1  I S  J 

OS  I  NO 

9  1  1  0 

3  1  30 

L  1ST 

USING 

9120 

3140 

LIST 

USING 

9130 

3 1  SO 

LIST 

USING 

9  1  AO 

3  1  AO 

1  1ST 

OS  1  NO 

9  1  SO 

3170 

1  1ST 

1 1S  T  NO 

9160 

Film  magnification  data  are  fit  to  a  specified  order  of 
polynomial  in  subroutine  LSQ.  The  coefficients  of  the  resulting 
polynomial  are  printed  and  stored  for  later  use. 


24 


4000  REM  -  READ  FIl M  MAGNIF  1C A r ION  DATA  AND  FIT  TO  A  2ND- ORDER  POLYNOMIAL 
4010  READ  Pi, Ml 
4020  FOR  1=1  TO  PI 
4030  READ  P<T>'Q<I> 

4040  NEXT  I 
4050  GOSUB  10000 
4060  M3  =  A ( 1 ? 1 > 

4070  M4-‘A<  2  »  1 ) 

4080  M5=A(3,1) 

4090  LIST  USING  9 1 70 , M3 »  M4 , M5 , E 


In  statements  4450  and  4480 ,  measured  values  taken  from  the 
film  (Figure  8)  are  used  to  compute  "instantaneous”  framing 
rates,  F3  and  F4 ,  for  successive  time  intervals.  These  computa¬ 
tions  are  tailored  to  16  mm  film  (40  frames  per  foot)  with 
timing  marks  applied  to  the  film  at  1  msec  intervals.  Successive 
"instantaneous"  framing  rates  are  averaged  to  provide  the  framing 
rate  data  which  are  fit  to  a  polynomial  in  subroutine  LSQ. 
Coefficients  of  the  polynomial  are  printed  and  stored  for  later 
use. 


4400  REh  •  READ  FILM  SPEED  DATA  AND  FIT  TO  A  2ND-  ORDER  POLYNOMIAL 
4410  READ  P1fM1»U<0) 

4420  FOR  1=1  TO  PI 
4430  READ  R(I)>U<I) 

4440  NEXT  I 

4450  F3=  ( <  II ( 2 )  ~U ( 1  )  )  *  1 0000 )  /  ( 3#  ( R  ( 2 )  R  ( 1 )  )  ) 

4460  FOR  1  =  1  TO  PI  -  2 

4470  P(I>  =  (lJ(Ifl)  U< 0 )  ) #10/3 

4480  F4=<  <U< I f2>  Ut Til >  >*10000>/<3*<ft< J+2)-R<l  +  l ) >  > 

4490  Q(I)s(F3+F4>/2 
4500  F3=F4 
4510  NEXT  I 
4520  P lap 1-2 
4530  GOSUB  10000 
4540  F5=A < 1 r 1 ) 

4550  F  6=A (2ft  ) 

4560  F7»A(3vl> 

4570  LIST  USING  9180»F5p F6>F7* E 


4800  REM  -  READ ,  COMPUTE  *  AND  PRINT  REAL  VALUE  OF  POINT  OF  IMPACT 
4810  READ  XO 

4820  X0= ( M3iM4*X0iM5*X0**2 ) *25  *  4 


4830  LIST  USING  9670 

4840  IF  SGN(XO)X)  THEN  4870 

4850  LIST  USING  9 1 90  *  ADS < XO > 


4860  GOTO  4880 


4870  LIST 
4800  LIST 
4890  LIST 
4900  LIST 
4910  LIST 


USING  9200 # ADS T XO) 
USING  9260 
USING  9270 
USING  9280 
USING  9290 


4920  LIST  USING  9680 
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Raw  data  for  the  pre-event  frame  are  processed  on  a 
section-by-section,  point-by-point  basis*  After  each  fringe 
number-measured  value  pair  has  been  read,  the  measured  value 
is  processed  to  compensate  for  film  magnification  errors  and 
to  convert  units  of  measurement.  Using  the  fringe  number  and 
X-coordinate  as  access  controls,  a  displacement  value  is  ob¬ 
tained  from  a  linear  interpolation  of  tabulated  fringe  coor¬ 
dinates.  These  initial  values  of  X  and  D  are  stored  and  then 
processed  in  subroutine  COR  to  determine  the  magnitude  of 
correction,  if  any,  which  must  be  applied  to  the  raw  values. 
After  the  corrections  for  film  magnification  have  been  made, 
the  raw  and  corrected  values  (and  a  count  of  the  number  of 
adjustments  made  in  subroutine  COR)  are  printed  at  the  terminal. 


‘ioou  id  n  id  mu  irini.i  nu.hii  hmiw  n  .*  n.uo  rok  h<e  ihpapi  pano  uma 

N0|0  I  OR  t  v  I  10  o  > 

sov'O  prinj 

PRUM  P  T  l  I.  V  I  I  0  > 

V'OIO  |  R  IO  I 

’  .<V  ,0  III  i J  I  ll:  .  Iil(;  94  10 
,:i  PK‘  J  0  T  1 1'  .*  INI.  V  I 
')<>  A)  PR*  I  NT 
“>000  RTOfi  I'l 
f>0°0  I  HI-  I  I  10  PI 

s  t oo  i\i  on  m  I  >  •r  u  - 

si  10  P <  i  »-  nivihUun  )  fhMK  i  » 

MV0 

fj  i  m )  ,M'  i  in  m  i 

S 1 40  t  OR  K  I  III  NVI  I  • 

mv>  u  xn.*.n  p<n  nn  n  mvo 
moo  Nr  XT  K 

S  I  VO  f  '0  •  (  X  <  R  I  y  J  >  P  l  I  >  )  /  <  X  (  R  I  *  I  >  X  (  K  y  J  >  ) 

'■>  i  no  o  ci  >  p  i  v  i  v  tpo*  <  n <  r  r  j ,  o  r  i  f  j  >  > 

K>  l  '  ■  0 

nvoo  omiiHi  tv, ooo 

MM  0  R  l  I  >  "  | 1  (  T  » 

?.>!■*  ‘0  Ut  1  y  I  V) -KM  I  ) 

fiV M)  Pin  MM  |  . 
f>v4o  m  :i  •  rv  >  r <  i  > 

“.vso  iKh=nc  n 
5VA0  HMMMKI) 

5V/0  PRINT  USING  9700, fM  1  )  ,&<!>, R<  l  )  ,S7,U<  I  )  y  YO 
f>2S0  NPXT  I 


The  corrected  X  and  D  data  pairs  for  each  pre-event 
section  are  fit  to  a  second  order  polynomial  in  subroutine 
LSQ.  Coefficients  of  the  resulting  polynomial  as  well  as  the 
measured  and  fitted  data  are  printed  at  the  terminal  and  the  user 


is  asked  if  the  next  higher  order  of  fit  is  desired.  A  positive 
response  results  in  a  repetition  (at  the  next  higher  order)  of 
the  curve  fitting  and  printing  process.  A  negative  response 
terminates  the  procedure  and  the  user  is  asked  to  specify  the 
order  of  fit  which  will  be  used  to  describe  the  particular 
section.  Pre-event  displacement  data  are  fit  to  a  polynomial  to 
provide  a  continuous  functional  relationship  which  is  used  during 
determinations  of  the  relative  displacement  of  the  panel  for  each 
of  the  varying  X  values  from  the  post  event  data. 


5400  R  EM  ~  FIT  X » D  MATA  TO  LEAST  SQUARES  POLYNOMIALS  ANM  PRINT  RESULTS 

5410  Ml =2 

5420  N4=0 

5430  GO SUM  10000 

5440  TE  N4--'t  THEN  5710 

5450  PRINT 

5460  PRINT 

5470  PRINT  USING  9430 f M t , p* i h 1  > 

5480  PRINT  USING  9220 

5490  FOR  J-0  TO  Ml 

5500  PR  TNT  US I NG  9230  .  J 9  A i Ji t  *  l  > 

5510  NEXT  J 

5520  PRINT  USING  9240*F 
5530  PRINT 

5540  PRINT  USING  9440 
5550  PRINT  USING  9450 
5560  PRINT 
5570  FOR  J=1  TO  PI 
5580  S0=A< 1 9 1 ) 

5590  FOR  K  "  1  TO  Ml 

5600  SO»SO  K  > 

5610  NEXT  K 

5620  PRINT  USING  *4AO.F<J>fRi  I )  0.1  <  ,J  i  •  50  9 1 1  f  J )  SO 
5630  NEXT  J 
5640  PRINT 
5650  PRINT 

5660  PR INI  “FIT  TO  NEXT  HIGHER  ORDER  POl  VNOMIAI  ?  iYLS/NO)"* 

5670  INPUT  X* 

5680  IF  X**“N0"  THEN  5710 
5690  M  l  -Ml II 
5700  GOTO  5430 
5710  PRINT 

5720  PRINT  “  DESIRED  ORDER  OF  FIT  FOR  USE  IN  EVALUATING  RFI  ATTUE  " 

5730  PRINT  “PANEL  DISPLACEMENT  FOR  POST  IMPACT  FRAMES  “ ; 

5740  INPUT  Ml 
5750  PRINT 
5760  PRINT 
5/70  GOSUD  10000 

5780  REM  STORE  COEFFICIENTS  OF  CUR ME  FIT  FOR  EACH  SECTION  LINE 
5790  FOR  1*0  TO  Ml 
5800  0< If 1 f 1 2  )  ~  A  ( I  f  1  • 1 > 

5810  NEXT  I 

5820  LIST  USING  9300* 12.0 •  I , 1 2 > . 0 < 2 9 12 ) 9 0 < 3 *  12 ) 9 0 < 4 9  1 2 )  90(59  12 > .E 
5830  0 < 10.12) ~M 1 
5840  0  <  9  9  1 2 ) WP 1 
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5850  Nl  a  I  I  :> 

58*0  RfM  coufi  i  rr  print  on i  of  i  ni  l  i vn.i 
5870  I  T  ST  US  I  NO  SO .  01  K  t  II  *  .  Mii  I  <  W  » 
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In  the  next  segment  of  the  program,  fringe  number- 
measured  value  pairs  for  all  post  event  frames  are  processed. 
Page  number,  frame  number,  shot  number,  and  the  time  after  the 
event  and  other  page  heading  information  are  printed  first. 
Processed  data  for  successive  sections  of  each  frame  are  then 
listed.  As  for  the  pre-event  case,  raw  measured  values  are 
corrected  for  film  magnification  before  they  are  listed.  Only 
the  corrected  values  are  listed,  however.  In  addition  to 
fringe  number  and  X  and  D  values  (in  device  coordinates)  the 
listing  includes  X  values  referenced  to  the  point  of  the 
event  and  displacement  values  referenced  to  the  pre-event 
position  of  the  target. 

7000  REM  READ  AMR  PRINT  AIL  FRAMES  01  RATA  hI'TER  PRF  IMPACT  FRAME 
7010  FOR  11-1  TO  NO  1 
7020  RF AO  N2 
70,10  PR  TNT  N.»; 

70 <10  I  1ST  USING  97*10 » »  P7 » NOT  I 
7050  F0=<N?/<F5Tf  A*Nr*TF7*N2''?>  **1000 
707.0  I  ..1ST  USING  77*20.  N?-T0 
7070  L  IST  USING  97*30 
7080  IT  ST  USING  97.40 
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7 0 9 0  I  1ST  US  I Nii  9A50 
7100  LIST  1 1ST  NS  9AA0 
7110  FOR  1 17“  l  rn  N< 

7120  REAM  PI 

7130  if  Pi  "O  I  MEN  7330 

7140  LIST  USING  9710>P* < 12+10) 

7150  f  OR  J>1  TO  PI 

71  AO  RE' All  ROMP'D 
M 7c  F  < I >  0(  |  ) 

7  1 00  P  <  I  )  i  fH  f  Hi  ■♦Pel  M  M5*P  (  1  >  2  >♦75  ,  4 
7170  J  T  l  IP  DM 
7200  FOR  K -  l  TO  N9t  I 

72  to  IF  X(K»J»  P<D  THEN  7730 
7220  NEXT  K 

7330  f ■*<>==  <  X  <  K  D  J  )  P  <  i  >  >  /  <  X  <  K  I  *  J  >  X  <  K  r  J  >  > 

7240  UU)-D(h  1  f  -J)+PO*<n<Kv  J>  IKK  If  PM 
7250  50 SI  IP  15000 
7 2 AO  S0~0  <  1  *  1 7  > 

72  70  FOR  I  =  1  TO  0<  JO* F7> 

7280  SO -SO  HI  <  I  -M  • 1 3  >  *  <  P  < i >  **t  > 

7290  NEXT  I 

7300  LIST  USING  9  A  90  *  F  <  T  >  *P«  I  >  *Q<  I  >  tPi  I  >~XO*G<  I  >--S0 
7310  NEXT  1 
7330  GOTO  7340 

7330  LTST  USING  9720 * P* < l 7 M 0 > 

7340  NEXT  12 
7350  P/  "P7U 

73  AO  Rf  M  GU  TO  NF  *  l  Mr  Aril  ANN  PAGE  OF  PRINTOUT 

7370  NEXT  II 

8000  PRINT  *  JOB  COMF1  I.  if  D  GOING  TO  E  ND  STATEMENT" 

BO  I  0  GOTO  20000 


Following  is  a  listing  of  the  various  image  statements 
used  in  conjunction  with  the  various  PRINT  USING  and  LIST  USING 
commands.  Content  of  these  statements  can  be  changed  to  suit 
the  heed  of  the  user.  Image  statements  listed  and  used  for  the 
sample  problem  were  used  for  an  aircraft  transparency  deflection 
measurement  application. 


9000 

9010 

9020 

9030 

9040 

9050 

90A0 

9070 

9000 

9090 

9100 

9110 


REM 

:i 

:o 

:o 


FOL 

OF 


IMAGE  STATEMENTS  FOR  TITLE  PAGE  PRINTOUT 

F 1 A-A  BIRDSTR IKE  PROGRAM 

###**♦*•**♦♦♦#♦*****♦##♦*♦#**#****# 
#»»#*##*»«***#♦##************ 
#*»#***####♦##♦•#*«**#♦*•*«♦#*****###***#### 
*♦*#****#♦♦*♦*##*♦**###**#************ 

LBS 

♦##  KNOTS 

OF  FILM  (AFTER  IMPACT)  EVALUATED  FOR  THIS  SHOT 


SHOT  NUMBER 
DATE  FIRED 
TRANSPARENCY 
PROJECTILE 
WEIGHT 
VELOCITY 
M  FRAMES 


LOWING  ARE 

the  form: 

Y 


THE  COEFFICIENTS  OF  LEAST- SQUARES  POLYNOMIALS 


»  AO  +  A1*X  + 


+  A N*X**N> 


29 


VI 20 
9130 
9  MO 
9150 
91  AO 
9170 

91  BO 
9190 
9200 
9210 
9220 
9230 

92  40 
9250 
9260 
9270 
9200 
9290 


:  FIT  TO  THE.  EXPERIMENTAL 

:  DATA.  FILM  MAGNIFICATION?  IN  REAL  INCHES  t  IS  GIVEN  AS  A  FUNCTION 
J  OF  MEASURED  INCHES  AND  FILM  SPEED r  IN  PICTURES  PER  SECOND t  13 
J  GIVEN  AS  A  FUNCTION  OF  FRAME  NUMBER  * 

SOCOEFFICIENTS:  AO  A1  A2  STD  ER  OF  ESI 

:  OF  TIM  MAGNIF:  *.*####'  !  !  !  #*#####'  '  •  1  #.#####!  !  !  '  ****** >  !  !  » 

SOFILM  SPEED!  #.#####»'»'  -#♦####*»!»'  *.*****<!!>  *.#*##*!!'» 

:  FUS  STA  140  LOCATED  ***.**  MM  10  THE  RIGHT  OF  CENTER  OF  MOIRE  DEVICE 
J  FUS  STA  140  LOCATED  ******  MM  TO  THE  LEFT  OF  CENTER  OF  MOIRE  DEVICE 
:  INSUFFICIENT  DATA  FOR  FIT  TO  A  *  ## -ORDER  POLYNOMIAL. 

J  PANEL  DATA  (REFERENCED  TO  THE  MlHRE  DEVICE  > . 

:  ac*;  -  • *  i* it ♦ 1 » ■  ' 

:  SID.  ERROR  OF  ESI  IMA  II  *.!****!  »  '  ' 

JO  PROGRAM  RUN  AT  ******  HRS  ON  ****** 

10  **  RESULTS  OF  CURVE  FITS  JO  PRE  ■  IMPACT  RANI  I.  SECT  IONS 

l  (REFERENCED  TO  IMF  MOIRE  DEVICE  > 

I OSEC  COEFFICIENTS  OF  LEAST  SQUARE  POLYNOMIALS  STD  ERROR 

:  NO  A  (  0  >  A  (  1  >  A  <,  2  >  At.  >>  AC  4;  OF  EST 


9300 
9400 
94  1  0 
9420 
9430 
94  40 
?  450 
9460 
9470 

9600 
96  I  0 
9690 
9630 
9640 
9650 
9660 
96  70 
9600 
9690 
9700 
9710 
9720 


*  #.***' i ' »  ♦,***' ! '  *  *.***«  ' !  » 

CEM  IMAGE  STATEMENTS  FUR  PRINTOUT  01  PRt 
FR  1  NGE  MEASURED  X  y  MM  D I SJ  lACF  MET  NT  >  MM 

NUMBER  RAW  l  ORR  RAW  CORR 

COEFFICIENTS  OF  A  *  ##  ORDER  POLYNOMIAL  FIT  TO 
I  RINSE  MEASURED  ACTUAL  CAL  CIJLATI  D 

NUMBER  X  f  MM  DTSPL.  r  MM  DISPL  » MM 

#*  **.***''ii  #*.#**#Mii 


*♦#**! ' ' '  *.#**'  «  >  ' 
IMPACT  RANI  l  DATA 

NO  TIMES  IN 
SLIER  COR 
THE  PRE  IMPACT 
DIFFERENCE 
ACT  CAT.  C 

it.****'  i  '  ' 


-■#♦*#!  !  '  ' 


: OFRAME  NUMBER  PRE  IMPACT 

REM  IMAGE  STATEMENTS  FOR  ROUTINE  DATA  PAGES 
:  I  SHOT  NUMBF.R  I*********************  ************ 

X OFRAME  NUMBER  **  <**.**#♦  MSEC  AFTER  IMPACT) 

:  X  COOR’D  I NA  Vf  *  DISPLACEMENT  «•  X  COORD  T NAIF  * 


F  RINGE  IN  MM  * 


I  N  MM  f 


I  N  MM  y 


NUMBER  W/R  TO  MOIRE 
TH  TRICE 


W/R  TO  MOIRl 
DEVICE 


W/R  TO 
FUS  STA  140 


PAGE  *#  OF  **  PACE 

DISPLACEMENT t 

IN  MM  r 

W/R  TO  PRESHOT 
POSI I  ION 


**  *******  *#♦*#*  ******  I*.*** 

#*  *11.**  **♦.**  ******  #*#.**  ** 
()*****#**###*###**####***#***#*******************#*#***##*+**#*#*# 
ODATA  NOT  AVAILABLE  IOR  ******  I*#**#****# I***#*##*#**#**##*******##** 


In  the  following  subroutine,  LSQ,  pairs  of  data  are  fit 
to  a  polynomial  using  a  least-squares  curve  fitting  routine. 

The  sets  of  and  data  are  used  to  form  the  following  series 
of  equations: 
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+  A1EXi 


+  •  • 


+  A  EX? 
m  1 


EXiYi 


A  EX. 

O  1 


Vxi 


+  A  EX?+1 
m  1 


EX?Y. 
1  1 


A  EX? 
o  l 


A. EX?+1 

1  l 


+  A  EX . 
m  i 


In  these  equations,  n  is  the  number  of  data  pairs,  m  is  the 
order  of  the  desired  fit,  and  A0 ,  Aj,  A2  ...  are  coefficients 
of  various  terms  in  the  equations.  Since  these  equations  are 
of  the  general  form. 


y  =  A„  +  Aj  X  +  A2X  +  .  .  . 

they  may  be  solved  simultaneously  to  yield  the  various  coeffi¬ 
cients.  In  subroutine  LSQ,  this  is  accomplished  by  forming 
the  various  summations,  placing  the  appropriate  values  in  arrays 
and  using  matrix  algebra  to  solve  for  the  coefficients.  The 
various  EX?  terms  are  formed  in  statement  10140  and  later 
stored  in  array  Y;  the  various  EX?Y^  terms  are  formed  in  state¬ 
ment  10180  and  stored  in  array  T.  Solution  of  the  coefficients 
is  accomplished  in  statements  10350  and  10360.  Upon  completion 
of  a  determination  of  the  standard  error  of  estimate  for  the 


fit  using  the  following  relationship. 


Std  Error  of  Estimate  = 


EY2  -  EA  x'V 
m 


control  is  returned  to  the  main  program. 


10000  REM  -  LEAST  SQUARES  SUBROUTINE 

10010  REM  -  UP  TO  50  BATA  POINTS  MAY  BE  PIT  TO  A  9TH  ORDER  POLYNOMIAL 

10020  M2=M1+1 

10030  IF  M2>P1  THEN  10430 

10040  MAT  A-ZER ( M2  f 1 ) 

10050  MAT  T~ZER<M2*1> 

10060  MAT  U=ZER(M2*M2) 

10070  MAT  Y~ZER<M2f M?) 

10080  FOR  J=0  TO  7*Mt 


10090  G  ( *J  )  “0 

10100  NEXT  J 

10110  S -Y2::s0 

10120  FOR  N  "  1  TO  F  I 

101 50  FOR  K= 1  TO  2* Ml 

1 0 1 40  G<  K  >  "-G  < K  1  if  <  N  >  ~K 

10150  NEXT  K 

10160  FUR  l-l  TO  M2 

10170  IF  L~:  t  THEN  10200 

10180  t  < I  f  I  )  ”T  (  L  f  .1  >  IO(NU  <  P (  N  > '  ( l  -  J  >  ) 

10190  GOTO  10210 

10200  TO  .  I  >  ::  I  ( I  y  I  )  HKN> 

10. MO  NFXT  t 

10770  Y2  Y7tU(N>''2 

10230  fi(0)-N 

10740  NEXT  N 

10750  I  Ml 

10  760  f  (JR  |  "  1  Ml  M2 

I  0  7  70  1=1.11 

10280  J~:  I  1 

10290  FOR  Kl  10  I 

10700  M-K  I 

10  UO  N  I  J 

10  >7"  t  Ml  •  0<M> 

10720  Ml  X  l  K 
10740  Ml  T  I 
10050  MAI  d  I NV  <  r  > 

10  160  MAI  A  :d*  I 
10370  FUR  II  ID  M3 
107,80  £>“::8  I  <  A(  T»M  KI  I  *  I  i  > 

10390  NEXT  1 
10400  Z~ABS<Y2  • .) ) 

10410  F>8UR'  Z/G<0  •  ) 

10420  C»U  I  D  10450 

10430  PRINT  US' TNG  92  I  0  *  M  I  «  I  t  <  ill  > 

10440  Ni  l 
10450  RETURN 


Film  magnification  correction  terms  were  determined  using 
data  taken  from  the  pre-event  frame-  When  the  target  is  de¬ 
flected,  a  change  in  the  magnification  of  its  image  occurs . 

Since  the  dark  band  used  to  define  a  measurement  point  is  formed 
on  the  surface  of  the  target,  evaluation  of  the  film  magnifica¬ 
tion  polynomial  using  post  event  measurements  may  yield  results 
which  are  somewhat  in  error-  Displacement  values  obtained  with 
use  of  these  measured  values  will  also  be  in  error-  The  mag¬ 
nitude  of  the  error  increases  as  the  post  event  measurement 
plane  is  displaced  in  either  direction  from  its  pre-event 
position- 

The  error  just  described  is  reduced  to  a  specified  mini¬ 
mum  value  in  subroutine  COR  using  an  iterative  correction 
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procedure.  A  graphical  illustration  of  one  step  of  the  cor¬ 
rection  procedure  is  presented  in  Figure  11.  In  this  procedure, 
use  is  made  of  the  fact  that  each  measured  value  lies  on  a  line 
which  makes  a  unique  angle  to  the  initial  position  of  the  target 
(or  the  X-axis).  With  respect  to  the  illustration  in  Figure  11, 
angle  a  should  equal  angle  B  for  this  condition  to  be  fulfilled. 
Angle  a  (A5,  statement  15020)  is  the  true  angle  which  the  line 
of  sight  makes  with  the  pre-event  target  surface.  Angle  6  (A6, 


Figure  11.  Illustration  of  Procedure  Used  to  "Correct"  Measured 
Values  for  Magnification  Errors. 


statement  15030)  is  the  angle  which  an  apparent  line  of  sight 
makes  with  the  pre-event  target  surface.  These  angles  are 
determined  in  subroutine  COR  using  device  constants  provided 
at  the  beginning  of  the  program  and  the  particular  X^  and  DA 
values  obtained  from  the  initial  interpolation  of  the  tabulated 
X  and  D  coordinates.  The  difference  between  these  angles  is 
compared  with  a  specified  minimum  acceptable  error  (statement 
15040).  If  the  difference  is  less  than  the  acceptable  error, 
control  is  returned  to  the  main  program. 

When  a  correction  must  be  applied,  the  magnitude  of  the 
correction,  AX  (S2) ,  is  determined  by  making  use  of  geometric 
relationships  illustrated  at  the  upper  left  of  Figure  11  and 
expressed  in  statement  15050.  The  correction  is  applied  to  the 
old  value  of  and  the  new  X^  is  used  to  determine  a  new  value 
of  by  interpolation  in  the  table  of  fringe  coordinates. 

Note  that  use  is  made  of  the  fact  that  the  meashred  point  must 
lie  in  the  center  of  the  dark  band  specified  as  the  fringe 
number  in  the  original  data  pair.  The  new  values  of  X  and  D 
are  used  to  compute  a  new  angle  8.  The  difference  between  a 
and  the  new  8  is  compared  with  the  minimum  acceptable  difference 
and,  if  necessary,  the  correction  process  is  repeated.  The 
examination  and  correction  procedure  is  repeated  until  the 
difference  between  these  angles  is  less  than  the  minimum 
acceptable  value. 

15000  ROI  SUBROUTINE  FOR  CORRECTION  OF  ERROR  DDE  TO  MAGNIFCATION 
15010  Y0~0 

1 5020  05--- AT N  (Rl/<  R3-P  <  I  )  )  ) 

15 030  A,  ATN  < <  R 1  Q (!))/( R3-P ( I ) ) ) 

15040  IF  ABS<A5-A6><5E~04  THEN  15140 

15050  S2*< ( <R3-P< I ) )/C0S(A6) ) *TAN < A5  A6 ) ) *SIN< A6 > 

15060  P<I)«P<I>+S2 

15070  FOR  K=*  1  TO  N9+1 

15080  IF  X(KfJXP(I)  THEN  15100 

15090  NEXT  K 

15100  PO-CXCK-l 9 J)-P<I ) 1 f J)~X(Kf J) ) 

15110  Q(  f  )=D(K-  1  *  J  >  1 1'O  t  (  D  <  K  f  J )  0  <  K  - 1  ?J>) 

15120  Y0-Y0+1 
15130  GOTO  15030 
15140  RETURN 


20000  END 


SECTION  V 

TEST  PROGRAM  AND  SAMPLE  PRINTOUT 


Except  for  order-of-fit  data  which  are  entered  from  the 
time-sharing  terminal  during  a  run,  all  data  used  during  the 
execution  of  PANDIS  are  read  from  DATA  statements.  Statement 
numbers  1  through  999  are  available  for  use  as  DATA  statements 
Following  is  a  listing  of  sample  data  which  are  read  from  these 
DATA  statements.  The  data  shown  in  these  statements  are  used 
to  produce  the  sample  printout  presented  in  this  section. 


An  explanation  of  the  data,  by  variable  name,  follows  this 
listing. 

l <>  (!.-> ( a  i vo  c..jv.>. iv*. r  '^n  ■<<? 

70  DA  1  A  I  •  »> I  0  ?  '*  ’<*  . 

SO  DAI  A  SI  f  ND  f  R  0  •  1 II »  1 1 1 »  (IN  III*  HI*  Iff 

31  DATA  |  X  |Ji:i  |  AI:MVI  Irll '^i,I  I  Til  rH 

37  DAT  A  SI  i  I  [Mil  I HUi Uhill  lr1IA»‘l  I  II  1 1 1 1 

33  DATA  snniON  1  INUI  HH  MW  I  Hi  ’At;  I  IViDI 

-10  DATA  ns  I  livCIDU  ri.iH  lODrD  4  l 1  *!?•> .  i  TICVI  D  LliU.S 

SO  DATA  l MAMMA Ui  DAI  I 

60  DATA  8  9  2 

/0  DA  1  A  1  * •  3  * 1  •  v  ,  f.J  ’  *  7  *  3  *  .  3  r  I  »  3  *  ♦  |  ♦  .  ’  ■ »  0  *  *  .'*• 

80  DATA  *  4  f  1  .  3  r  .  *73  f  7  ,  5  .  *  873  .  3  ,  3 

90  DATA  I  4  *  7  f  4  ♦  4  4 

100  DATA  0*0.  I  .  .86*2*  I  .77.  ^ .  7  4  *  >>  .  4  4  *  3 *  4  .  3f6f  5.  1  7 

1 10  DATA  7r6.0Jr8r6.9-9F7.77f i Of  8. 64- I  I  . 1 ?f 10 *  >7 .13.11*  3  1 
130  DATA  » 895 
140  DATA  7 

145  DATA  3f1.6f  2f1.185f  -  1  f  »  715f  0"  1  f  Of  .255*  If  .66f -2*  *85 

ISO  DATA  7 

133  DATA  3r  l  .5f~2p  i  *005»~  I  •  .  6.1  b  f  0  f  0  •(>  f  .355r  1  •  .76f~2f-  .95 
160  DATA  7 

165  DATA  -3f  1  .4f- 2f  .?85f  1  f  .515f0f  •••• .  I  f  0  f  -  .  455  f  -  J  f  -  *  86  f  -2  f  I  .  OS 

170  DATA  1 

ISO  DATA  S 

1 85  DATA  -3 f 1  . 6  f -  2  f 1 ♦ 29  f  1 f . 92 fOf . 55  f 1 f . 01 f 1 f  ♦ 23 

190  DATA  Of  .65 r~ l*  *83 
200  DATA  8 

205  DATA  -3 f 1 . 5 f -2 f! .19f  1 f . 87  f  0  f . 45  f 1 f  ~  . 09  f 1 f  • . 35 
210  DATA  Of  . 75 f~ If -.93 
220  DATA  8 

225  DATA  ~3f  1 .4f -2f  1  f09f-1  f  .72f0f  .35f  1  »-♦  19f  1  f~.45 

230  DATA  0f-.85f~1f~1.03 

240  DATA  2 

250  DATA  0 

260  DATA  9 

265  DATA  ~3f 1 . 59f ~2f 1 . 29f  -1 f . 9Bf0f .68f 1 f . 37f2f - . 02 
270  DATA  2  f - ♦ 57  f 1 f - ♦ 79  f  0 f  ~ ♦ 93 
280  DATA  9 

285  DATA  3 f 1 . 49 f -2f 1 . 19f ~1 f . 88 f 0 f . 58f 1 f ♦ 27f 2f ~ ♦ 12 
290  DATA  2  f  •  67  f  1  f  ~  •  89  f  0 f  —  1  •  03 


r 

i 


Statement 

Number 


Variable  Name/Description 


Rl,  R2,  N9 ,  R8,  R9 
Al,  A2 ,  FI,  F2 ,  N3 

Statements  10  and  20  contain  the  six  moire  device 
parameters  and  other  values  used  to  establish 
limits  for  subsequent  computations  and  procedures. 
Note  that  N3  has  a  value  of  3,  indicating  that 
profiles  are  to  be  obtained  along  three  sections 
of  each  frame. 

P$(l)  to  P$ ( 9 ) 

P$(ll) 

P$  (12) 

P$ (13) 

Statements  32  through  36  are  descriptions  of  the 
section  lines  which  were  used  for  each  frame. 

S$,  D$ ,  P$(0) 

P$  (10)  ,  Y8 ,  Y9 ,  NO 

Statements  40  and  50  contain  descriptions  of  test 
target  conditions  and  a  specification  of  the 
number  of  frames  which  will  be  analyzed  for  the 
particular  test. 

PI,  Ml 

The  number  of  points  and  the  order  of  fit  desired 
for  the  film  magnification  correction  polynomial. 

P(l) ,  Q ( 1 ) ,  P  (2) ,  R(2) , . . . 

P  (6 )  ,  Q  ( 6 )  ,  P  ( 7 )  ,  Q  (7)  ,  .  .  . 

Data  for  film  magnification  determination 

PI,  Ml,  U(0) 

The  number  of  points,  the  order  of  fit,  and  the 
distance  to  the  first  frame  of  the  event.  These 
values  are  used  with  the  film  speed  data  to 
determine  the  coefficients  of  a  polynomial  which 
describes  the  framing  rate  of  the  camera. 

R (1 ) ,  U(l) ,  R ( 2) ,  U  ( 2) ,  .  .  . 

R(8) ,  U (8) ,  R(9) ,  U (9 ) , . . . 

Data  for  film  speed  determinations 


Statement 

Number 


130 


140 


145 


150,160 


155,165 


170,250 


Variable  Name/Description 


X0 

Distance  from  origin  of  moire  device  coordinate 
system  to  point  of  interest  on  target 

PI 

Number  of  data  pairs  for  first  section  in  pre¬ 
event  frame 

U(l) ,  R(l) ,  U ( 2 ) ,  R(2) , . . . 

Fringe  number,  measured  value  pairs  for  first 
section  in  pre-event  frame 

PI 

Number  of  data  pairs  for  remaining  sections  in 
pre-event  frames 

U (1) ,  R(l) ,  U  (2) ,  R(2) , . . . 

Fringe  number,  measured  value  pairs  for  section 
in  pre-event  frame 

N2 

Frame  numbers 


180,200,220 
250,260,280  PI 

Number  of  data  pairs  for  each  section  in  the 
various  post  event  frames.  Note  that  a  value  of 
zero  was  entered  in  statement  250  indicating  data 
were  not  available  for  this  section. 


All  other 
statement 

numbers  Q(l),  P(l),  Q(2),  P(2),... 

Data  pairs  for  sections  of  all  post  event  frames. 


The  next  three  pages  are  printout  that  was  received  at  the 
the  terminal  during  the  run  using  the  sample  data.  Only  results 
of  the  analysis  of  certain  pre-event  data  are  printed  at  the 
terminal.  These  data  are  printed  to  permit  the  user  to  select 
best  fit  data  for  use  during  analysis  of  subsequent  post  event 
frames.  In  each  of  the  pages  the  user  was  asked  for  a  response 
after  results  of  the  first  curve  fit  were  printed.  As  shown. 
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the  response — YES — resulted  in  another  printout  and  question. 

The  negative  response  resulted  in  a  request  for  specification  of 
the  order  of  fit  to  be  used  with  the  data  that  were  presented 
at  the  beginning  of  the  printout.  A  similar  sequence  was 
followed  for  the  additional  sections  of  the  pre-event  frame. 

At  the  end  of  the  third  page,  the  number  of  frames  which  had 
been  processed  and  a  notification  that  the  job  had  been  com¬ 
pleted  are  printed  to  signify  that  work  on  the  run  has  been 
terminated. 


SAMPLE  TERMINAL  PRINTOUT 


Pirst  Section 


♦RUN 

>  USED  AS  ARRAY  &  SCAL  AR  NAME  5 E » H » M »N » W 


SECTION  1-INCH  ABOVE  IMPACT  POINT 


FRINGE 

MEASURED 

X  f  MM 

DISPLACEMENT  *MM 

NO  TIMES  IN 

NUMBER 

RAW 

CORR 

RAW 

CORR 

SUBR  COR 

-3 

101*88 

100.86 

-24.02 

-25.41 

1 

82*30 

79 . 80 

--14.58 

1 7 . 05 

2 

-1 

53.76 

50.98 

-8.80 

-9.65 

2 

0 

6*21 

6.84 

1.14 

1*26 

2 

0 

-26*51 

23.47 

4.88 

4.32 

3 

1 

68.54 

68.37 

0.35 

0*26 

1 

....  *? 

-89.99 

-93.41 

-6  ♦  72 

3*23 

1 

COEFFICIENTS  OF  A  2  ND  ORDER  POLYNOMIAL  FIT  TO  THE  PRE  IMPACT 
PANEL  DATA  (REFERENCED  TO  THE  MOIRE  DEVICE). 

A ( 0  >  *  J.&5370E+00 

A ( 1 >  =  -1 .03775E-01 
AO)  =  -1.&7575E-03 
STD.  ERROR  Or  ESTIMATE  *  8.3R858E-01 


FRINGE 

MEASURED 

ACTUAL 

CALCULATED 

DIFFERENCE 

NUMBER 

X  *  MM 

DISPLvMM 

DISPLvMM 

ACT  -  CAI  C 

-3 

10 ♦ 08 6 Li 01 

-25 . 4073E+00 

25 *859 IE TOO 

45. 1782E  02 

-2 

79 . 8 OOF:  >00 

-1 7*0539Ft00 

1 7 . 3286ET00 

27.4673E  02 

1 

50.983E+00 

-9A.4530E* 01 

-79. 929 IE -01 

-16.5230E  01 

0 

6S.370E  01 

12.5744E  01 

86.5055E- 02 

39.1503E-0? 

0 

23.472ET00 

43. 1690F  01 

31.66301:  01 

11 .506 IE -01 

1 

•  68 ♦ 372E f 00 

25.5676E  02 

91 •5397F- 02 

-65  *9  7200  02 

....  2 

93.414L 100 

•  37 .31 70 E  01 

-32.7503E-01 

43. 32A4F-03 

FIT  TO  NEXT  HIGHER  ORDER  POLYNOMIAL?  <YES/MO)?YES 


COEFFICIENTS  OF  A  3  RD  ORDER  POLYNOM IAL  Ell  TO  THE  PRE  IMPACT 
PANEL  DATA  (REFERENCED  TO  THE  MOIRE  DEVICE). 

A ( 0 )  =  1 . 73528E+00 

A ( 1 )  *  -1.22948E-01 
A<2)  =  I . 70283E-03 
ACT)  =  2.58884E  OA 

STD.  ERROR  OF  ESTIMATE  =  7.O7723E-0I 


FRINGE 

MEASURED 

ACTUAL 

CALCULATED 

DIFFERENCE 

NUMBER 

X*  MM 

DISPL *MM 

DISPLvMM 

ACT  CALC 

-3 

1 0  *  086ET0 1 

-25 « 4073E+00 

-25 • 3306ET00 

-76  *  72 1 2E -03 

-2 

79  *  880E1 00 

-17 ♦ 0539E  f 00 

-17.A317E TOO 

57 • 7805E-02 

-  1 

50 ♦ 983ET00 

96.4530E  01 

86. 161  IE  01 

10.2919E  01 

0 

68 . 370E-  01 

12  * 5744L  01 

81  .591 5E -02 

44. 1524E  02 

0 

-  23 • 472E  f 00 

43.1690E  01 

36.4950E  01 

66.7 404E-02 

1 

-68 • 372E+00 

25.5676E02 

13.5380E  Ot 

10*  981 3E -01 

....  n 

-93  »  41 4ET 00 

-32*  3170E-01 

-37* 4904E-01 

51.7338E  02 

FIT  TO  NEXT  HIGHER  ORDER  POLYNOMIAL?  <YES/NO)?NO 

DESIRED  ORDER  OF  FIT  FOR  USE  IN  EVALUATING  RELATIVE 
PANEL  DISPLACEMENT  FOR  POST  IMPACT  FRAMES  ?3 
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■>V\T  V^*-.  ±3*'* 


SAMPLE  TERMINAL  PRINTOUT  -  Second  Section 


SECTION  THROUGH  IMPACT  POINT 


FRINGE 

MEASURED 

X.MM 

DISPLACEMENT »MN 

NO  TIMES  IN 

NUMBER 

RAW 

CORK 

RAW 

CORR 

SUBR  COR 

-3 

97.64 

95*55 

■29*77 

32.61 

—2 

76.79 

72  *  57 

-20.20 

■•24 . 35 

3 

-1 

46.82 

42  *  82 

-10*91 

12.12 

3 

0 

•  2 . 62 

-2.40 

0.48 

0.44 

1 

0 

-36.42 

32*02 

6.70 

5.89 

3 

“1 

-79.70 

■75.68 

7 . 1 6 

4.08 

3 

....  o 

101*73 

99.00 

5.26 

2.47 

2 

COEFFICIENTS  OF  A  : 

>  ND  ORDER 

POLYNOMIAL 

FIT  TO 

THE  PRE  IMPACT 

PANEL  DATA  (REFERENCED  TO  THE  MOIRE  DEVICE ) ♦ 

A (  0  )  «  8.45184C  02 

Ad)  *  -1.92726E  01 
A (2)  *  -1.68482E--03 
STD*  ERROR  OF  ESTIMATE  a  9.99634E-01 


FRINGE 

MEASURED 

ACTUAL 

CALCULATED 

DIFFERENCE 

NUMBER 

X  f  MM 

DISPLfMM 

DISPLfMM 

ACT  -  CALC 

•3 

95 . 546ET00 

-32. 609 5E +00 

-33.8796E+00 

12* 7014E- 01 

72. 56 6E TOO 

-24 ♦ 3455ET00 

22  *  9420E+00 

■•  14.0352E  01 

•  l 

42 • 825ET00 

- 12*  U84ETOO 

11  •  4279F+00 

■••■69 «  049 IE- 02 

0 

-24.050E > 0 1 

44 ♦ 2338E  02 

36*9241 E  02 

73 ♦ 0975F-03 

0 

“32 • 01 6ET00 

58 ♦ B826E  01 

43.5886E-0 1 

15 . 294 OE  01 

1 

75 ♦ 676E LOO 

40.8426E  01 

48.5154E  01 

76.7280E  02 

....  “) 

98. 998E+00 

24.713 4E -01 

24.8279F-  01 

U.4441E-03 

FIT  TO  NEXT  HIGHER  ORDER  POI  YNOMIAL?  (YES/NO)? YES 


COEFFICIENTS  Of  A  3  RD  ORDER  POLYNOMIAL  FIT  TO  THE  PRE  IMPACT 
PANEL  DATA  ( REFERENCED  TO  THE  MOIRE  DEVICE), 

A < 0 )  ®  - 1.351 81 E™ 01 
A<1)  =  -2.20498E-01 
A <2)  =  --1.66065E-03 
A (3)  a  3.73797E-06 
STD.  ERROR  OF  ESTIMATE  =  7.54936E-01 


FRINGE  MEASURED 
NUMBER  X  f MM 


ACTUAL 

DISPLfMM 


CALCULATED  D 1 1  TERENCE 

DISPL.MM  ACT  CALC 


3 


95  *  546E TOO  -32  * 6095E+00 


33 « 1Q26E+00 


49 * 3179E-02 


72  *  566E+00  24 . 3455E+00 


-23  *  4523E+00  -89 . 32 1 9E-02 


1 

0 


42 * 825E+00  12.11 84E  TOO 

-24 . 050E-01  44 . 2338E-02 


12.3299ET00 

38*5457E~02 


2 1 • 1589E-02 
56*  881 7E -03 


0  -32.016ET00 


58  * 8826E-01  50 . 9945E- *01 


78  *  8801E-02 


1  --75.676E+00 

2  98*  990E+00 


40#  8426E-01 
24  # 71 34E-01 


54 • 2099E-01  -13.3673E-01 

1 7  #  9 1 69E-0 1  67 . 9655E -  02 


FIT  TO  NEXT  HIGHER  ORDER  POLYNOMIAL?  (YES/NO >TNO 

DESIRED  ORDER  OF  FIT  FOR  USE  IN  EVA!  HATING  RELATIVE 
PANEL  DISPI  ACEMENT  FOR  POST  IMPACT  FRAMES  *3 
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SAMPLE  TERMINAL  PRINTOUT  -  Third  Section 


SECTION  1-INCH  BELOW  IMPACT  POINT 


FRINGE 

MEASURED.  X>MM 

DISPLACEMENT. MM 

NO  TIMES  IN 

NUMBER 

RAW 

CORR 

RAW 

CORR 

SUBR  COR 

-3 

93*10 

89*54 

-35.92 

-40.75 

2 

-  2 

70*98 

65*42 

-25*06 

-27 ♦ 56 

n 

-  1 

39.57 

34*17 

13*11 

14.74 

3 

0 

-11.75 

-10*57 

2.16 

1*94 

2 

0 

-46.64 

•40*57 

9*17 

7*46 

4 

1 

91  ♦  15 

-81.51 

19*61 

9.13 

2 

....  o 

-113.77 

104*38 

24*48 

7.97 

3 

COEFFICIENTS  OF  A 

2  ND  ORDER 

POLYNOMIAL 

FIT  ro 

THE  PRE  IMPACT 

PANEL  DATA  (REFERENCED  TO  THE  MOIRE  DEVICE > ♦ 

A  <  0 )  «  -1 ♦ 56133E+ 00 
A ( 1 )  *  -2 «  81 086E  0 1 
A< 2 )  a  1 « 80357E-03 
STD  *  ERROR  OF  ESTIMATE  =  6.90966E 01 


FRINGE 

MEASURED 

ACTUAL 

CALCULATED 

DIFFERENCE 

NUMBER 

X  ?  MM 

DISPL t MM 

DISPL f MM 

ACT  -  CALC 

-3 

89 ♦ 540E+ 00 

-40.7533F f 00 

41 . 1898EI00 

43  *  6462E-02 

65  *  422F TOO 

27 • 5633LM00 

27 ♦ 6699E TOO 

1.0  ♦  661 3E-02 

i 

34-1  69E  f 00 

14*  7423E  L00 

1 3 «  271 3E TOO 

•14 ♦ 7100E-01 

0 

-10*  566ET00 

19 * 432 IE  01 

12. 0 7 21 E  01 

73-599 7E  02 

0 

•40  *  567E+00 

74.6073E  0.1 

68.7330E-01 

58.7427E  02 

-1 

-81 . 507E+00 

9J ♦ 265  4  P  01 

93 « 6728E- 01 

-24.0739E-02 

....  '■» 

~i0.43&E+0i 

7?.??7?r:  01 

91.2773E  01 

.1 5  *  4938E-02 

FIT  TO  NEXT  HIGHER  ORDER  POLYNOMIAL?  ( YES/NO >?YES 


COEFFICIENTS  OF  A  3  RD -ORDER  POLYNOMIAL  FIT  TO  THE  PRE  IMPACT 
PANE!  DATA  (REFERENCED  TO  THE  MOIRE  DEVICE ) ♦ 

A  <  0  )  *-  1  »  71 183E  LOO 

Ad)  «  -2 * 99229E  01 
A ( 2 )  a  :*1  .  7433 IE -03 
A ( 3 )  «  2.52445E  06 

STD*  ERROR  OF  ESTIMATE  *  5.28980E-01 


FRINGE 

MEASURED 

ACTUAL 

CALCULATED 

DIFFERENCE 

NUMBER 

X .  MM 

DISPL. MM 

DISPL t MM 

ACT  -  CALC 

•3 

89 ♦ 540 E  LOO 

-40.7533F400 

40. 6694EF00 

-83.8928E  03 

65 . 422F100 

77  *  563 3E f 00 

- 28.O425F+00 

47.9263E  02 

1 

34. 169E+00 

-  14 -7423F+00 

-  1 3  *  8706E \ 00 

■:?7. 1632F -02 

0 

-10.546E+00 

1  9, 4321 E  0.1 

12. 521 6E  01 

69, 1051T- 02 

0 

-40.567E+00 

74.6073E  -01 

73.8944E  “01 

7 1.289 IE-03 

X 

-81  *  507E+00 

91.2654E-01 

97 ♦ 2894E-01 

-60.2397E-0? 

....*7 

-10 . 438ET01 

79 ♦ 7279E- 01 

76. 5639E-01 

31.6407E-02 

FIT  TO  NEXT  HIGHER  ORDER  POLYNOMIAL?  <YES/NO)?NO 

DESIRED  ORDER  OF  FIT  FOR  USE  IN  EVALUATING  RELATIVE 
PANEL  DISPLACEMENT  FOR  POST  IMPACT  FRAMES  ?3 


1  2  JOB  COMPLETED-  GOING  TO  'END'  STATEMENT 

* 
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The  next  four  pages  present  the  sample  printout  produced 
at  the  high-speed  printer.  Presentation  of  these  pages  completes 
instructions  for  the  use  of  the  moire  device  and  program  PANDIS. 


rifc-A  BI^OSTPIKe  PROGRAM 


SHOT  NUKbC* 


TEST  PRObLEM 


DATE  FjPLC 


OCTObEP  4  1980 


TRANSPARENCY 


CURVlD  LINES 


project  ill 

*C ICHT 
VELOCITY 


IMAGINARY  u  A  L  L 
C.uC  LfS 
3  kNGTS 


2  FRAMES  of  FILM  liriEP  IMPACT)  LwAlUATLC  fcr  this  shot 


FOLLOWING  APE  THE  COEFFICIENTS  OF  LEAST -SO J  A WE  S  POLYNOMIALS 
uF  THE  FORM: 

Y  =  AO  ♦  A  1  *  X  ♦  ...  ♦  /  N  *  Y  *  *  \  f 

FIT  TO  THE  EXPERIMENTAL 


,ATA.  FILM  vAGNIFTCA  T I  ON  f  IN  °lAL  INCHES,  is  GIVEN  AS  A  FUNCTION 
OF  MEASURED  INCHES  AN  D  FILM  Sr  E r  D  ♦  IN  PICTURES  PER  SLCONOt  IS 
oiVEN  AS  A  FUNCTION  OF  FRAME  MJ M  b  E  R  • 


COEFFICIENT: 

A  C 

A  1 

A2 

STO  CR  Of  ES 

F I L  v  M A  ON  IF  : 

-i  .o2°7SE*LJl 

3.c25l?i>GH 

T.O^Zf'OL-O? 

FIL"  SFLED: 

?.  i>S027G*{j  J 

1  .?4S54C«UD 

-2.78473t-02 

6.  928?0f  *00 

FUS  STA  1  4  J  LOCATED  *5.24  TO  THE  n IGHT  OF  CFNTER  OF  MOIRE  DEVICE 

**  PLSULTS  OF  CURVE  FITS  To  PPE~lMPACT  PANEL  SECTIONS  ♦♦ 
(REFERENCED  TO  THE  %TRl  DEVICE* 

iEC  COEFFICIENTS  OF  L E A S T -S wU A »■ t  POLYNOMIALS  STD  ERROR 


NO 

A  (0* 

All) 

A  12) 

A  1  3  ) 

A  (  4  ) 

OF  EST 

1 

1.  7  2  5  E  ♦  T  C 

-1  *2i9F.-Cl 

-  1  • 7T  3r -0  3 

2*539E-06 

0*  COCE ♦QO 

7.08E-D1 

c. 

-1. 3S2E-C1 

-2 .2D5**-C1 

-1 .t61E-n3 

3.726E-C6 

0 • COOE ♦OU 

7.S5E-31 

S 

-1. 7 1 2E ♦ DO 

-2.9S2E-G1 

• 1 •  74  3E-C  3 

2.524£-r6 

0«  C03E ♦  03 

5.29E-01 

PROGRAM  RUN  AT  I6U0S  hRS  ON  llOSFu. 


SAMPLE  PRINTOUT  -  Title  Page 
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SHOT  NUMBER 


TEST  PROBLEM 


PAGE  2  OF  t  PAGES 


FRA  HE  MJHPEP  --  PRE-JTPACT 


X-C00R3INATE 

,  DISPLACEMENT  , 

X-CoOPDlNATE , 

DISPLACEMENT , 

r  RINSE 

IN 

IN  *M, 

1'4  MMf 

IN  M*% 

NUMBER 

J/P  TO  KOIkE 

a /R  TO  M0  IRE 

m/e  TO 

*/R  TO  PPE SHOT 

DEVICE 

DEVICE 

f US  STA  140 

POSITION 

SECTION 

1-INCH  AtOVE 

IMPACT  POINT 

-  r 

1CJ. 653 

-  2  ^  •  4  J  7 

196.095 

-3.077 

_  -N 

79.88j 

-17.054 

175.117 

0.576 

-  1 

50.963 

-9.64« 

146.220 

-1.029 

n 

6.637 

1.257 

102.074 

0.442 

-21.172 

4.317 

71.765 

J.667 

-1 

-6c. 272 

0.256 

26.865 

-1.09b 

.  'I 

-93.414 

-2.232 

1.823 

u.517 

SECTION 

THRCLGh  Impact  point 

-  t 

95.5^6 

-32.610 

1 9C • 7b  3 

3.493 

-? 

72.666 

-24.746 

lb  7 .8  03 

-C.893 

-1 

42.625 

-12. IIP 

1 3  6  •  r6  2 

C  .  2  1  2 

£ 

-2.4C5 

r.442 

92.622 

C  •  C  r  7 

-32.016 

5.638 

63.221 

u.789 

-1 

- 75 • t  76 

4  .C  64 

19.561 

-1.327 

-2’ 

-96.998 

2.471 

-3.761 

0.660 

SECTION 

1-INCH  6 tLC* 

IMPACT  POINT 

-  ”5 

8  9  •  5  4  C 

-40.75? 

164.777 

-0.06.4 

—  7 

65.422 

-27.5b? 

lb  0 • 6  59 

D.479 

-1 

34.169 

-14.742 

129.405 

-C.B72 

- 1  Li  .  566 

1.943 

o4.67l 

C.691 

t 

-45.567 

7.461 

54.670 

0  •  07 1 

-] 

-el. 507 

9.127 

13.73C 

-0.6C2 

-  2 

-104 .364 

7.97? 

-9.147 

3.316 

SAMPLE  PRINTOUT  Cont'd  -  Pre-event  Frame 


SHOT  NUMBER  --  TEST  PROBLEM  PAGE  3  OF  4  PAGES 


FRAME  NUMBER  --  1  ( 

C . 34  70 

MSEC  AFTEP 

IMPACT  ) 

X-CGORDZNATE 

,  DISPLACEMENT, 

X -COORDINATE  , 

DISPLACEMENT  , 

FRINGE 

IN  MM  * 

IN  MM, 

IN  MMt 

IN  M  M  f 

NLMtER 

W/P  TO  MCIRL 

W  /  ft 

TO  MOIRE 

«/P  TO 

b/R  TO  PR f  S HO T 

DEVICE 

DEVICE 

F  US  STA  140 

POSITION 

oECTIOI. 

1-INCH  AbOVt 

IMPACT 

POINT 

ICC.fcSo 

-2?  .4j? 

196. T95 

-C.077 

_  -> 

66.57b 

-1C  .214 

lbl  .615 

9.779 

-1 

66.c63 

-P.661 

lb?. 120 

12.469 

0 

4 5 • 006 

.27^ 

140.243 

15.269 

1 

1 u  .  52 1 

21 .f 37 

1U5.75F 

20.761 

1 

-12  .  164 

19.462 

o  3  •  Oti  3 

1  6  •  5  r  9 

'• 

-54.031 

1 3.9U0 

40. f C6 

10.953 

“1 

-79. *65 

7.441 

15.26? 

c.Ges 

SECTION. 

through  impact  point 

-3 

9S.S«fc 

-32.6  IP 

19C.76  3 

^.453 

-? 

PG.213 

-16 .709 

175. 455 

5.671 

-  1 

5V.ll 1 

-6  .  c67 

ii>4.?47 

11.17b 

37.3C3 

6  .*6P 

132. r4P 

17.33c 

1 

<..016 

4.0. 317 

97.251 

20.9C3 

1 

- 1  V  .  b  1'j 

21.437 

75.627 

1 7.84^ 

■) 

-61 .u9C 

1^.722 

34.147 

11.437 

-  1 

-65.0*1 

1  3  .P2P 

IP. 146 

c.  720 

oECTiCN 

1-INCH  cELO* 

I  mpact 

PC  1M 

_  ? 

8  5  .  0  4  J 

-4P.757 

164 .777 

-U.0C-4 

72.  V? 3 

-  c  4 .14° 

16B  .  170 

7  •  6  p  1 

- 1 

51.372 

r27 

146. 609 

11.61b 

2  9  .  2  5  7 

6  •  3  b  1 

12  4 .494 

1 7.27o 

i 

-  b  .  S  1  7 

lr  .r6l 

b6  .  72C 

1  5  *  2  0  7 

i 

-27. L£  1 

27. 7 2 P 

o?  .  1  66 

! o.271 

-67.61b 

21  •  r  2  4 

27.722 

11.757 

-I 

-  0  L-  •  U  1  U 

1° . 7o4 

5.227 

9. 1P7 

SAMPLE  PRINTOUT  Cont’d  -  Frame  #1 


SHOT  NUMBER 


TEST  PROBLEM 


PAGE  **  OF  <4  PAGES 


FRAME  NUMBER  --  7  <  C.6RJ8  M<EC  AFTER  IMPACT) 


X-COORDINAU  f 

DI SPL  A  CFHTNT  , 

X-COQPOINATE, 

01 6PL ACEh£NT  v 

FRINGE 

I*  HH  • 

IN  « 

IN  MHt 

IN  hHf 

NLMPER 

ta/P  TC  HOIRL 

W /k  TC  HOTRf 

m/P  TO 

m/P  TO  PRESHOT 

DEVICE 

DEVICE 

fus  sta  mo 

POSITION 

CATA  NfT 

AVAILABLE  Fup 

SECTION  1-TnCH 

ABOVE  IMPACT  PCINT 

SECTION 

through  impact 

PCINT 

-  :■ 

100.37? 

-  2  6  •  T  6  6 

195.609 

9.  mi 

86.578 

-lf.rm 

It  1. PIS 

19.033 

-1 

71 .oP9 

2.9  17 

166.336 

25.281 

r* 

55.706 

19.589 

15  1  •  r  3  3 

31.596 

l 

39.623 

25 .9iT 

139.760 

37.197 

t 

1  7 . 5  C  9 

it  .c71 

112.796 

91.C56 

J" 

-25.2*6 

37.927 

69 .9449 

33*609 

1 

-52.o99 

2°. 965 

*2.293 

23.667 

r 

-72.5P7 

29.525 

22.650 

It. 635 

section 

1-INCH  PtLC.  IMPACT  POINT 

»  A 

99.9*1 

-33.375 

190.218 

1  C  .  3  2  2 

-  7 

8  C  .  2  1c 

-16 .709 

175.965 

1 6.922 

-1 

63.610 

159.097 

23.997 

W 

9  7  •  c  1  9 

9.75] 

1 H  2  •  p  56 

29.3*2 

1 

31.192 

2  2  •  7 1  7 

126.928 

36.962 

“ 

1  G  •  L  9  9 

36. c  73 

10c  .286 

91.963 

? 

-31.91b 

9n  .  Ct'' 

62.721 

39 .C°9 

1 

-5c. 518 

33  .A -7 

36  .719 

2  9. 36  G 

7 

-7-5.771 

2  0  •  1  c  e 

16.966 

16.377 

SAMPLE  PRINTOUT  Concluded  -  Frame  #2 
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a  Moire  Frinqe  Technique, 


AFWAL-TR-81-3005 ,  Air  Force  Wright  Aeronautical  Laboratories, 
Wright-Patterson  Air  Force  Base,  OH,  May  1981. 
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DATE 

ILMED 


